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Abstract 

In this paper, the general plane problem for an infinite strip con- 
taining multiple cracks perpendicular to its boundaries is considered. 
The problem is reduced to a system of singular integral equations. Two 
specific problems of practical interest are then studied in detail. 

The first is the investigation of the interaction effect of multiple 
edge cracks in a plate or beam under tension or bending. The second 
problem is that of a rectangular plate containing an arbitrarily 
oriented crack in the plane of symmetry. Particular emphasis is placed 
on studying the problem of a plate containing an edge crack and subjec- 
ted to concentrated forces. The plate has the dimensions of a standard 
compact tension specimen and is intended to simulate CTS. 


1 . Introduction 

The plane problem for an infinite strip containing an edge crack 
simulating a single edge notch specimen, a beam, or a plate and the 
problem of a rectangular block with an edge or an internal crack simula- 
ting a compact tension specimen are two of the more widely studied 
geometries in fracture mechanics. Aside from the standard finite ele- 
ment methods (e.g. [1]), a wide variety of analytical and numerical 
methods have been used to solve the problem. Some of the significant 
techniques used in these studies are the Wiener-Hopf method (e.g. [2]), 
the method of weight functions (e.g. [3]), the method of Laurent series 
(e.g. [4,5]), the conformal mapping technique (e.g. [6]), the method 
of integral equations (e.g. [7,8]), and the method of boundary 


*This work was supported by NASA-Langley Research Center under the Grant 
NGR 39-007-011 and by NSF under the Grant ENG 78-09737. 
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collocation (e.g. [9-11]). In this paper the basic problem of multiple 
cracks for an infinite strip is considered by using the method of singu- 
lar integral equations. The paper has two objectives. The first is to 
provide the solution of the interaction problem for a beam or a plate 
containing two or three edge cracks perpendicular to the boundary and 
subjected to membrane loading or pure bending. The second is to give 
an analytical solution to the crack problem for a rectangular block or 
a compact tension specimen which is subjected to arbitrary crack surface 
tractions or concentrated (body) forces. 

2. Integral Equations of the Problem. 

The basic crack geometry under consideration is shown in Figure 1. 
It is assumed that y=0 is a plane of symmetry with respect to loading 
as well as crack geometry and the conditions of the plane strain or the 
generalized plane stress are satisfied. In addition to arbitrarily dis- 
tributed crack surface tractions, the medium may be acted upon by 
arbitrarily located concentrated forces P shown in the figure. First, 
it may easily be shown that(*) for a pair of point dislocations with 
densities g and h located at the point (x 0 ,y 0 ) and defined by 

[v(x,y o +0) - v(x,y o -0)] = g(x 0 ,y 0 )<$(x-x 0 ) , 

Cu(x,y o +0) - u(x,y o -0)] = h(x 0 ,y 0 )6(x-x 0 ) (la,b) 

the stress state in an infinite plane may be expressed as 

a xx ( x,y ) = 9( x o»yo) G xx( x >y» x o*yo) + h(x 0 ,y 0 )H xx (x,y,x 0 ,y 0 )> 

° yy (x,y) = g(x Q ,y 0 )G yy (x,y,x 0 ,y 0 ) + h(x 0 ,y 0 )H yy (x,y ,x 0 ,y 0 ) , 

o xy (x,y) = g(x 0 ,y 0 )G xy (x,y,x 0 ,y 0 ) + h(x 0 ,y 0 )H xy (x,y ,x 0 ,y 0 ) , (2a-c) 

where u and v are the x and y-components of the displacement and the 
influence functions are given by 

v 'For example, by using complex potentials [12] or standard Fourier 
transforms. 


* 


G xx (x,y,x 0 ,y 0 ) 
H xx (x,y,x 0 ,y Q ) 
G yy (x,y,x 0 ,y Q ) 
H y y(x,y,x 0> y 0 ) 
G xy (x,y,x 0 ,y 0 ) 
H X y (x,y,x 0 ,y 0 ) 


_ 2u 

(x 0 -x)[(x 0 -x) 2 - (y-y 0 ) 2 ] 

tt(1+ic) 

[(x„-x) 2 + (y-y 0 ) 2 ] 2 

_ 2y 

(y-y 0 )C(y-y 0 ) 2 + 3(x 0 -x) 2 ] 

it ( 1 +k ) 

[(x 0 -x) 2 + (y-y 0 ) 2 ] 2 

_ 2 y 

(x„-x)[3(y-y 0 ) 2 + (x 0 -x) 2 ] 

IT ( 1 +IC ) 

C(x 0 -x) 2 + (y-y 0 ) 2 ] 2 

2y 

(y-y 0 )[(y-y 0 ) 2 - (x 0 -x) 2 ] 

tt(1+ic) 

[(x 0 -x) 2 + (y-y 0 ) 2 ] 2 

2y 

(y-y 0 )C(y-y 0 ) 2 - (x 0 -x) 2 ] 

TT ( 1 +K ) 

C(x 0 -x) 2 + (y-y 0 ) 2 ] 2 

_ 2u 

(x 0 -x)[(x 0 -x) 2 - (y-y 0 ) 2 ] 

ir(l+<) 

C(x 0 -x) 2 + (y-y 0 ) 2 ] 2 


(3a-f) 


In (3) p is the shear modulus and k= 3-4v for plane strain and < = 
(3-v)/(l+v) for generalized plane stress, v being the Poisson's ratio. 

Similarly, for a pair of concentrated body forces P and -P (per 
unit thickness) acting in y direction at points (m,n) and (m,-n), respec- 
tively, the stresses in the plane may be expressed as follows: 


<J XX (x,y) = P Q xx (x,y,m,n), 

a yy (x,y) = P Q yy (x,y,m,n), 

oJ y (x,y) = P Q xy (x,y,m,n), (4a-c) 


where 

Q xx (x ’ y ’ m ’ n) = 2V('l+<) { (x-mJ^T (y-n) 2 Ck_1 

y+n r . 4(x-m) 2 t, 

(x-m) z + (y+ n) z ” (x-m) 2 + (y+n) 2 ^ 
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Let us now consider the stress state in an infinite strip 0<x<H 
parallel to the y axis for which y=0 is a plane of symmetry. Using 
Fourier transforms and the conditions of symmetry it may be shown that 

4(x,y) S~ {[a(A,+xA 2 ) + ^ A 2 ]e'“ x 

0 

+ [a{A 3 +xAi,) - — g -- Aiijs 0 *} COSay da, 

o* (x.y) = r {MAj+xA,) + ^ A 2 ]e -ax 

JJ Q 

+ [a(A3+xA4) - K ^ A4]e a } COSay da, 

a^y (x,y) = - ^ /" {[a(Ai+xA 2 ) + A 2^ e aX 

-[a(A 3 +xA 4 ) - ~ 2 ~~ A 4 ]e } sinay da, (6a-c) 

where A 1 ,...,A 4 are unknown functions of a and are determined from the 
boundary conditions at x=0 and x=H. It is clear that if one imposes 
the boundary conditions 

a (o,y) = o (o ,y) = o (H,y) = ^ (H,y) = o, -~<y<“ , (7) 

XX Xy aa Ay 

on a solution for which the stress state is given by 

a io -(x,y) = a^.(x,y) + a!?j(x,y) + a^.(x,y), (i,j) = (x,y), (8) 
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where a !?. is the sum of the stresses due to symmetric dislocations located 

* J 

at points (x 0 ,y 0 ) and (x 0 ,-y Q ) as given by (2} with g(x 0 »y 0 ) = g(x 0 ,-y 0 ) 
and h(x 0 ,y 0 ) = -h(x 0 >-yo)» one v/ould obtain the formulation of a strip 
0<x<H, -«<y<» which is free of surface tractions, is subjected to con- 
centrated forces P and -P, and contains the dislocations defined by 
(1) at (x 0 ,y 0 ) and (x 0 ,-y 0 ). If the dislocations mentioned are the 
only defects in the strip, then by substituting from (2-6) and (8) into 
(7) and by inverting the Fourier transforms, one can determine the 
unknown functions A*,..., At* in terms of P, g, and h and thus, obtain 
the closed form solution of the problem. Needless to say, since the 
concentrated forces and the dislocations are "point functions", the 
solution for any number of symmetric forces and dislocations may.be 
obtained by a simple superposition of suitable solutions given by (2-8). 
Also note that the more general nonsymmetric solution may be obtained 
by expressing the strip solution (6) in terms of infinite Fourier trans- 
forms and by eliminating terms involving (y+n) in (5). 

Instead of dislocations if the strip contains cracks along c<x 0 = 
t<d, y 0 = B = constant under a given set of surface tractions, by inte- 
grating the solution found for the dislocations in t one would obtain 
a system of integral equations for the unknown density functions g and 
h which may now be considered as functions of t only. Referring now to 
Figure 1, let g(t), h(t), and g(t), -h(t) be the density functions 
defined by (1) for the cracks II and III, respectively for which c<t<d, 
y Q = + B. Also, let the strip have an additional crack I along a<x 0 = 
t<b, y 0 = 0. Considering the symmetry of the problem, for crack I one 
may write 

jj- [u(x,+0) - u(x,-0)] = h(x) = 0 , 

~ v(x,+0) v(x,-0) = f (x) , a<x<b. (9a, b) 

In addition to concentrated forces ^P let the’ strip be subjected to 
the following crack surface tractions 
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(10a, b) 


a yy (x, 0) = Pi(x), a xy (x,0) = 0, a<x<b , 
tf yy ( x » B ) = a yy (x,-B) = p 2 (x) , 

tf xy (x,B) = -cr xy ( x ’“ B ) = P 3 OO 1 c<x<d . (11a, b) 

Thus, the terms a*? . in the superimposed stress state given by (8) 

* J 

must contain the contributions from crack I with the density functions 
g=f and h=0, crack II with g and h, and crack III with g and -h. 

Following the procedure obtained above, after some straightforward 
manipulations the functions A-j(a), (i=l,...,4) may be determined from 
the boundary conditions [7] as follows: 

Ai(a) = (2aD) [4a 2 H 2 -(i<-l )Lj]R]i + [4a 2 H 2 -(<+l )1-2]R2 

+ [(1-k)L 3 + 2aH(e a ^*-ice 01 ^)]R 3 + [(1+k)L 3 - 2aH(e a ^+Ke a ^)]R 4 }, 

A 2 (a) = D" 1 [L 1 R 1 +L 2 R 2 +(2aHe aH +L 3 )R3 + (2aHe aH -L 3 )R 4 ] , 

Aa(a) = (2aD) \ [4a 2 H 2 +(<-l )L 4 ]Ri + [(k+ 1 )Ls-4a 2 H 2 ]R 2 

+[(k- 1 )L3+2aH(ice aH -e aH )]R 3 + [(1-Hc )L 3 -2aH(<e aH +e aH )]R 4 }, 

A 4 {ct) = D" 1 CL 4 R 1 +L 5 R 2 +(L3+2aHe" aH )R3 + (L 3 -2ctHe aH )R 4 ] , (12a-d) 

where 

D(a) = e 2ctH + e“ 2aH - 4a 2 H 2 - 2 , (13) 

and the functions L-j(a) ( i=l .... ,5) , and R-j(a), (i=l,...,4) are given 
in Appendix A. Thus, the complete solution of the problem is obtained 
once the density functions f, g, and h are determined. One may note 
that because of the assumed symmetry in formulating the problem the 
condition (10b) is automatically satisfied. Substituting then from 
(8) giving the combined stresses in the strip into the crack surface 
boundary conditions (10a), (11a) and (lib) one would obtain the system 
of integral equations to determine the functions f, g, and h as follows: 
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7 f* C 4 — + k n (x,t)]f(t)dt + 7 / k 12 (x,t)g(t)dt 
7 /k 13 (x,t)h(t)dt * [P k 14 (x) + p x (x)], a<x<b , 

| ^ k 2 i(x,t)f(t)dt + 7 /^t^x + k 22 (x»t)]g(t)dt 

+ 7 / k 23 (x,t)h(t)dt = -^f" CP k 24 (x) + p 2 (x)], c<x<d, 

7 / k 31 (x,t)f(t)dt + 7 / k 32 (x,t) g(t)dt 
7 T a 77 c 

+ ff / Ct=x + k 33(x,t)]h(t)dt = [Pk 34 (x) + p 3 (x)], C<X<d, 

d 

(14a-c) 

where the kernels k...(x,t) and the functions k^, (i,j=l,2,3) are given 
in Appendix B. 

If the cracks are internal cracks as shown in Figure 1, then the 
solution of the integral equations (14) must satisfy the following single- 
valuedness conditions: 
b 

/ f (t)dt = 0 , 
a 

d 

/ g(t)dt = 0 , 
c 

d 

/ h(t)dt = 0 . (15a-c) 

c 

For the internal cracks the integral equations (14) have ordinary 
Cauchy kernels and their solution may be obtained in a simple manner 
by using, for example, the technique described in [13]. After deter- 
mining the density functions f, g, and h, the stress intensity factors 
at the crack tips may be defined and evaluated as follows: 

Ma) - ^x)o yy (x, 0 ) = SZZ&n f(x), 

Mb) - >^Tx-BT Oyy (x, 0 ) = - ^ ^ f(x) , 
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MO - c yy (x,B) • g(x), 

mo = m a xy( x »B) = ^ ^ s£Gzr mo, 

MO = ’i3 ^dT c yy (x,B) = - |^ g(x) , 

k 2 (0 = o xy (x,B) = - h(x) - 06a-f) 

where k-| and ^ are Modes I and II stress intensity factors. 

For edge cracks, the asymptotic behavior of the kernels and the 
nature of the singularity of the solution of the integral equations were 
treated in [14], In the problem under consideration the asymptotic 
analysis of the kernels have been performed for the limiting cases of 
(a=0, b<H), (a>0, b=H) , (c=0, d<H), (c>Q, d=H), and (c=Q, d=H). 

Omitting the analytical details, only the results are given in equations 
(14) and Appendix B. From (14) and Appendix B it may be seen that the 
singular part of the kernels on the main diagonal of the system of 
integral equations, for example, for a=0 and c=0 is given by 

k s ( x »t) = t>x ' t+x + (t+x)' 1 " (t+x) d ' 

The kernel k g (x,t) is a generalized Cauchy kernel. The peculiarity of 
this particular generalized Cauchy kernel is that even though k $ (0,t)=0 
for 0<t<H and k g (x,0) = 0 for Ckx<H , the part of k s excluding the Cauchy 
kernel, i.e., the sum of the last three terms in (17), becomes unbounded 
as x ‘and t approach the end point zero simultaneously, that is, for t->0 
with x=At, A being an arbitrary nonzero constant. Because of this property 
it was shown in [14] that for a=0 f(a) and for c=0 g(c) and h(c) would 
be bounded. This is also the result one would expect on physical grounds. 
The technique for solving the edge crack has also been described in [14]. 

3. The Results and Discussion 
(a) The Infinite Strip. 

The first problem considered is the problem of interaction of two 
symmetric internal cracks in an infinite strip under uniform tension. 
Referring to Figure 1 and to the integral equations (.14)’, for this 
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crack geometry we have a=b, c=H-d, P=0, p 2 (x) = -a Q , and p 3 (x) = 0. 

Some calculated results are given in Table 1. The stress intensity 
factors obtained for 1 / H = 0.05, £=(d-c)/2 are indistinguishable from 
the values given in [15] which are obtained for an infinite plane. Note 
that as the crack distance 2B decreases ^ also decreases and k 2 becomes 
more significant. The angle e shown in the table is a measure of the 
probable crack growth direction in brittle materials and is obtained 
from the simple assumption that along this direction a 0Q (r, 0 ) is maxi- 
mum, where r « H-d [16]. Here 9>0 means that the cracks would propagate 
away from each other. 

Some results for the problem of interaction of two symmetric edge 
cracks in a strip under uniform tension or pure bending away from the 
crack region are given in Figures 2-5. For this geometry we have a=b, 
c=0, and d<H. The figures also show the values of kj for a single edge 
crack (as the dashed line). For a single edge crack k 2 and 0 are zero. 
Again note that is smaller than the corresponding single edge crack 
value, k 2 becomes more significant as B decreases, and in brittle mater- 
ials the cracks would tend to propagate away from each other. 

The results for three edge cracks in a strip under uniform tension 
or bending are shown in Figures 6-11. In all cases k 2 (d)<0 meaning 
that the outside cracks would again grow away from the crack in the 
middle. It may be observed that for a strip under tension and pure 
bending the Mode I stress intensity factor for the middle crack is less 
than that of the outer cracks which, in turn, is less than the corre- 
sponding single crack value. Comparing the results given in Figures 
2-5 with those given in Figures 6-9, it may be seen that the presence 
of the middle crack "relaxes" the stress state in the strip resulting 
in smaller stress intensity factors. Figure 10 shows the effect of the 
crack depth on the stress intensity factors. For b=d«(B,H) the results 
reduce to the value for a single edge crack in a semi- infinite plane, 
i.e., k x = 1.1216, k 2 = 0. As crack lengths increase, first the inter- 
action effect and then the free boundary effect (from x=H) dominate. 
Consequently the values of kx(b) and k x (d) somewhat decrease before 
increasing sharply. Figure 11 shows an example for a beam or plate 
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Table 1. Stress intensity factors in a strip containing 
two symmetric internal cracks, £=( d-c)/2. 


m 

B /l 

W °</1 

k 2 / a 0 /e 

9 (°) 


0.5 

0.7797 

- 0.1175 

16.430 


1.0 

0.8512 

- 0.0616 

8.194 


1.5 

0.9052 

- 0.0308 

3.887 

0.05 

2.0 

0.9395 

- 0.0163 

1.992 


5.0 

0.9953 

- 0.0001 

0.157 


10.0 

1.0053 

- 0.00001 

0.014 


20.0 

1.0060 

0.0000 

0.000 


0.5 

0.7992 

- 0.1199 

16.363 


1.0 

0.8749 

- 0.0624 

8.076 


1.5 

0.9310 

- 0.0307 

3.774 

0.1 

2.0 

0.9660 

- 0.0162 

1.920 


5.0 

1.0219 

- 0.0001 

0.106 


10.0 

1.0247 

- 0.00001 

0.003 


20.0 

1.0248 

0.0000 

0.000 


0.5 

0.8846 

- 0.2570 

15.578 


1.0 

0.9749 

- 0.0656 

7.634 

OJ 

• 

o 

1.5 

1.0437 

- 0.0330 

3.648 

2.0 

1.0839 

- 0 . 0 T 55 

1.641 


5.0 

1.1096 

- 0.0001 

0.019 


10.0 

1.1097 

0.0000 

0.000 
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under three point bending^*). Here the stress intensity factor for the 
outer crack becomes zero for B=4H (the load distance), decreases as B 
approaches zero and hence goes through a maximum for a certain value 
of B/H which may depend on the crack depth. 

(b) The Compact Tension Specimen. 

The problems of the rectangular plate and the compact tension speci- 
men are solved by letting c=0 and d=H in the basic strip problem (see 
Figure 1). In this case both ends of the outer cracks are treated as 
if they are the free ends of an edge crack. The integral equations 
(14) with the generalized Cauchy kernels are still valid. In all the 
problems considered for this geometry, it is assumed that the tractions 
p 2 and p 3 on the surfaces of the outer cracks are zero. 

The calculated stress intensity factors for a square plate con- 
taining an internal crack and subjected to uniform tension perpendicular 
to the crack are shown in Table 2. Here the parameters e and x repre- 
sent the eccentricity in the crack location and the relative crack 
length, e=0 being the symmetrically located central crack. A limited 
comparison of the stress intensity factors calculated in this paper and 
that obtained in [5] for the crack geometry used in [5] is given in 
Table 3. In [5] only ki(b), the greater of the two stress intensity 
factors, is given. The agreement appears to be quite good. 

The stress intensity factor in a square plate containing a 
syrnnetrically located central crack and subjected to concentrated wedge 
forces P (per unit thickness) in the middle of the crack is shown in 
Table 4. The table also shows the results given in [11]. Again, the 
agreement appears to be very good. Note that for the infinite plate 
the result is k x = P/tt/E where l = (b-a)/2. The effect of the location 
of the concentrated wedge force on the stress intensity factor is shown 
in Figure 12 for a specific crack geometry. 

7*1 I 

' ' In Figure 11 the strength of material solution is used for the 
uncracked strip. See [7] and [8] for the effect of concentrated 
forces and through thickness stress distribution obtained from the 
elasticity solution on the stress intensity factors. 
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Table 2. Stress intensity factors in a uniformly loaded square 
plate containing an internal crack. c=0, d=H, B=H/2, 
e=(b+a-H)/H (eccentricity), x=(b-a)/H, Pi(x)=-a , 
P2 = 0, P3=0 (see Fig. 1 and eq. 11) u 


e 

\ 

ki(a)/a Q /(b-a)/2 

ki(b)/CT 0 /(b“-a)/2 


0.1 

1.0140 

1.0140 


0.2 

1.0554 

1.0554 


0.3 

1.1233 

1.1233 


0.4 

1.2162 

1.2162 

0 

0.5 

1.3339 

1.3339 


0.6 

1.4810 

1.4810 


0.7 

1.6774 

1.6774 


0.8 

1.9914 

1.9914 


0.9 

2.7129 

2.7129 


0.1 

1.0139 

1.0140 


0.2 

1.0553 

1.0559 


0.3 

1.1234 

1.1269 

0.1 

0.4 

1.2172 

1.2326 


0.5 

1.3370 

1.3919 


0.6 

1.4880 

1.6669 


0.7 

1.6972 

2.3186 


0.1 

1.0151 

1.0161 


0.2 

1.0596 

1.0692 

0.2 

0.3 

1.1345' 

1.1759 


0.4 

1.2469 

1.3870 


0.5 

1.4278 

1.9113 


0.1 

1.0246 


0.3 

0.2 

1.0975 

1.1527 


0.3 

1.2431 

1.5170 

0.4 

0.1 

1.0789 

1.1227 
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Table 3. Comparison of the stress intensity factors calculated 
from (14) and that given in [5] for a rectangular plate 
under uniform tension. Pi=-a , a=0.545H, b=0.755H, 



£=(b-a)/2. 


V 



B/H 

ki(a)/a Q /e 
(Eq. 14) 

0.25 

1.2000 

0.50 

1.0617 

0.71 

1.0406 

1.0 

1.0382 

GO 

ki(b)/a 0 /£ 
(Eq. 14) 

1.2061 

1.0655 

1.0518 

1.0507 

1.0507 

ki(b)/a Q /e 

[5] 

1.205 

1.066 

1.050 

1.050 

1.0507 


Table 4. Stress 
central 
a=H-b, 

intensity factor 
crack which is 
£=(b-a)/2 , Pi (x) ; 

in a square plate containing a 
loaded by wedge forces. B=H/2, 
=-P5(x). 

£/H 

k ; /(P/*/£) 
(Eq. 14) 

k!/(P/W T) 
(Ref. [11]) 

0.05 

1.02907 

1.0279 

0.10 

1.11278 

1.1115 

0.15 

1.2514 

1.2499 

0.20 

1.4437 

1.4418 

0.25 

1.6889 

1.6866 

0.30 

1 .9921 

1.9894 

0.35 

2.3804 

2.3772 

0.40 

2.9561 

2.9523 

0.45 

4.1717 

4.1665 
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The results for a rectangular plate containing an edge crack and 
subjected to uniform tension are given in Figure 13 and Tables 5 and 6 
(see insert in Figure 13). The figure shows the effect of the crack 
length on the stress intensity factor for two values of B/H. The effect 
of length- to-width ratio B/H of the plate on the stress intensity fac- 
tor is given in Table 5. For B/H = 2 the result is practically the 
same as that found for an infinite strip with an edge crack. Table 6 
shows the effect of the crack length on the stress intensity factor in 
a square plate containing an edge crack and subjected to uniform ten- 
sion (see insert in Figure 13). The table also shows the results given 
in [6]. It is seen that the agreement between the two sets of values is 
very good. 

Table 5. Stress intensity factor in a uniformly stressed 
rectangular plate containing an edge crack. The 
crack length b=0.5H, a=0, Pj(x) = -c Q . 

B/H 0.25 0.5 1.0 2.0 

MbJ/oQ/F 4.8298 3.0103 2.8263 2.8254 2.8250 


Table 6. Stress intensity factor in a uniformly stressed 
square plate with an edge crack. a=0, B=0.5H, 
Pi (x) = -a Q . 


b/H 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

ki(b)/a 0 /ET 
(Eq. 14) 

1.2295 

1.4877 

1.8483 

2.3245 

3.0103 

4.1525 

6.4044 

12.0013 

ki(b)/a Q /b 
(Ref. [6]) 

1.23 

1.49 

1.85 

2.32 

3.01 

4.15 

6.40 

12.0 


The results for a square plate with an edge crack and subjected 
to concentrated wedge forces P (per unit thickness) on the crack sur- 
faces are given in Figure (14). Here the variable is the location of 
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the wedge force. As m-*b the stress intensity factor would tend to the 
infinite plate value given by 

ki(b) = — — — . (18) 

ir/b-m 

That is, it becomes unbounded. Also, as ithO due to the "bending" effects 
the stress intensity factor would again increase. Hence, k x (b) goes 
through a minimum for a certain value of the load distance m. Similar 
results have been observed in [8] for an infinite strip with an edge 
crack which is subjected to a concentrated wedge force. 

The remaining results in this paper concern an edge-cracked rec- 
tangular plate having the overall dimensions of a compact tension speci- 
men which is under a pair of concentrated (body) forces P only. The 
geometry of the CTS is given by Figure 15. Referring also to Figure 1, 
it may be seen that W=0.8H, B=0.48H, m=0.2H, and n=0.32H. First, one 
may note that the effect of n, the load distance in y-direction on the 
stress intensity factor is rather insignificant. This may be seen from 
Table 7. As seen from Table 8, this is not the case for the load dis- 
tance m in x-direction. It may be observed that for large values of m 
the stress intensity factor becomes negative. This is due to the fact 
that in the uncracked plate under concentrated loads +P the elasticity 
solution would give compressive stresses on part of the y=0 plane. 


Table 7. The effect of load distance n in y-direction on 
the stress intensity factor in a compact tension 
specimen under concentrated forces P. B=0.48H, 
m=0.2H, b=0.6H (see Figure 1). 


n/H 

ki(b)/(P/b/H) 


0.32 0.3 

7.8487 7.8489 


0.2 0.1 

7.8492 7.8479 


0.05 0.01 

7.8480 7.8257 


0 

7.8549 
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Table 8. The effect of load distance m in x-direction on the 
stress intensity factor in a compact tension speci- 
men under concentrated forces P. B=0.48H, n=0.32H, 
b=0.6H (Figure 1) 


m/H 0 0.2 

k 1 (b)/(P/E/H) 10.2859 7.8492 


0.4 0.6 0.8 

5.3738 2.7207 0.4255 


0.85 1.0 

-0.02137 -1.4287 


Table 9. Stress intensity factor for the compact tension 
specimen under concentrated forces P. m=0.2H, 
n=0.32H, B-0.48H, k*=k x (b)/(P/b/H) ( Figures 1 and 15). 


b/H 

0.50 

0.52 

0.54 

0.56 

k* 

6.1535 

6.4122 

6.7042 

7.0341 

b/H 

0.66 

0.68 

0.70 

0.72 

k* 

9.6294 

10.4387 

11.3997 

12.5507 


0.58 

0.60 

0.62 

0.64 

7.4130 

7.8487 

8.3525 

8.9385 

0.74 

0.76 

0.78 

0.80 

13.9435 

15.6493 

17.7694 

20.4519 


The stress intensity factor for the compact tension specimen 
corresponding roughly to the standard load location m=0.2H, n=0.32H 
(Figures 1 and 15) is given in Table 9. The same results calculated 
in terms of the standard A/W values shown in Figure 15 and normalized 
with respect to K 0 =P/v^W are given also in Table 10 (second column). 

For practical applications the stress intensity factors given in Table 
10 may be approximated by, for example, a polynomial through a least 
square curve fit in the following form: 


- £. / A ,H 
~ /W W 


N 

Z 

n=0 


(A/W)' 


(19) 


It should be noted that K which appears in (19) 
standard Mode I stress intensity factor used in 
is related to the stress intensity factor Mb) 


and in Table 10 is the 
fracture mechanics and 
defined in (16) by 


K = fit Mb) . 


( 20 ) 
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Table 10. 


Stress intensity factor for the compact tension 
specimen under concentrated forces P. n=0.4W, 
m=0.25W, B=0.6W, H=1.25W, K^P/vtf (Figures 1 and 15). 


A/W 

K(b)/Ko 
(Eq. 14) 

K(b)/K 0 
(Eq. 19) 

K(b)/K Q 
(Ref. 10) 

0.375 

6.8981 

6.9272 

6.8177 

0.400 

7.3304 

7.3182 

7.2787 

0.425 

7.8102 

7.7844 

7.7829 

0.450 

8.3449 

8.3264 

8.3396 

0.475 

8.9501 

8.9462 

8.9603 

0.500 

9.6381 

9.6499 

9.6591 

0.525 

10.4263 

10.4484 

10.2989 

0.550 

11.3364 

11.3598 

11.3643 

0.575 

12.4020 

12.4111 

12.2337 

0.600 

13.6465 

13.6427 

13.6541 

0.625 

15.1204 

15.1007 

15.1132 

0.650 

16.8832 

16.8594 

16.8569 

0.675 

19.0155 

19.0041 

18.9659 

0.700 

21.6283 

21.6442 

21.5518 

0.725 

24.8794 

24.9146 

24.7719 

0.750 

29.0000 

28.9787 

28.8558 


N 

Table 11. Coefficients 
c o C 1 c 2 

c n which 
c 3 

appear in 
c 4 

Equation (19). 
c 5 c 6 

3 

-43.252 

336.265 

-709.559 

529.552 




4 

171.451 

-1261.74 

3664.38 - 

■4693.42 

2298.41 



5 

-207.428 

2006.98 

-7369.26 

13530.5 

-12437.1 

4668.26 


6 

27.606 

-91.036 

44.302 

448.366 

-630.114 

-134.171 

512.814 

7 

23.663 

-70.262 

38.850 

310.498 

-406.739 

-45.894 

107.933 


- 17 - 


Table 11 gives the coefficients c n of the polynomial (19) for various 
values of N. For N=6 the (approximate) results calculated from (19) 
are given in Table 10 (column 3) along with the values calculated from 
the expression given in [10] (column 4). These results show that for 
a compact tension specimen, particularly in analyzing fatigue crack 
propagation data, the expression given by (19) or that developed in 
[10] may be used with a certain degree of confidence. 
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APPENDIX A 


The functions L-j and which appear in (12): 

Lj = -e 2ctH - 2aH + 1, L 2 = e 2aH - 2aH - 1, L 3 = e aH - e' aH 
L 4 = e" 2aH - 2aH - 1 , L 5 = e‘ 2ctH + 2 a H - 1 . 


R.j (a) = F.. (a) + (a) + H.j(a) + P.j(a), i = 1 

■, b f 

FlU) = - yyy / a f(t) at e dt , 

, d _ f 

G 3 (a) = - ^yy COSaB / g(t)at e dt , 

Hi (a) ‘= ~xt sinaB / h(t)e a ^(at+l)dt, 

K 1 c 

n / \ _ e" 0111 sinan (K-l-2am)' n 

p i {a) = jrtptI P ’ 

F 2 (o) = - ~y/f(t) e" at (l-at)dt, 
l d _ f 

G 2 (a) = - yyy COSaB / g(t)e~ at (l-at)dt, 

1 d t 

H 2 (a) = - yyy sinaB f h(t)at e~ at dt, 

n / \ _ e aITI sinan(-ic-l+2am) D 

2\ a l 4p ( k"+1 ) 5 


F 3(a) = - 


f(t)a(t-H)e" a(H " t) dt, 


G 3 ( a ) = 


H 3 (a) = 


P 3 (a) = 


- ^yy COSaB /g(t)a(t-K)e" a ^ H " t ^dt, 
-yr SinaB I h(t)e” a(H " t) [a(H-t)+l ]dt, 

Kfi c 

e sinanrK-l-2a(H-m)l p 

4i T^Tj p » 


(A. 1 ) 
(A. 2) 
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[l-a(H-t)]dt , 


F„(a) - - ^ 

G 4 (o) * - ^7j- cosaB /^g(t)e" a ^ H_t ^ [l-a(H-t)] dt, 

H 4 (o) = - sinaB / hUMt-Hje'^'^dt, 

p%(.) ■ ffia i ' c * 1 ~ 2 ' ,(l * ,)] p • 


(A. 3) 
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APPENDIX B 


The kernels k^-Cx.t) and the functions k i4 (x), (i,j=l,2,3): 

k u (x,t) = G(x,t) + G(x-H,t-H) + k fl (x,t) , 
k 12 (x,t) = 2CQi(B,x,t) + Q 2 (B,x,t) - Q 2 (B,H-x,H-t) + k gi (x,t)] 

k 13 (x,t) = 2[Q 3 (B,x,t) + Q 4 (B,x,t) + Q 4 (B,H-x,H-t) + k hl (x,t)] , 

k 2 i(x,t) = Qi(B,x,t) + Q 2 (B,x,t) - Q 2 (B,H-x,H-t) + k f2 (x,t) , 

k 22 (x,t) = G(x,t) + G(x-H.t-H) + Qi(2B,x,t) + Q 2 (2B,x,t) 

- Q 2 (2B,H-x,H-t) + k (x,t), 

k 23 (x,t) = Q 3 (2B,x,t) + Q 4 (2B,x,t) + Q 4 (2B,x,t) + k h2 (x,t) , 
k 31 (x,t) = Q s (B,x,t) + Q 6 (B,x,t) + Q s (B,H-x,H-t) + k f3 (x,t) , 
k 32 (x,t) = Q 5 (2B,x,t) + Q 6 (2B,x,t) + Q 6 (2B,H-x,H-t) + k g3 (x,t) , 

k 33 (x,t) = G(x,t) + G(x-H,t-H) + Q 7 (B,x.t) + Q 0 (B,x,t) 

- Q 8 (B,H-x,H-t) + k h3 (x,t) , (B.l) 

where, 

G(x,t) - - ^ . (B-2) 


Qi(B.x.t) = • 

Qa(B,x,t) - ’ 

Q 3 (B.x,t) = ’ 
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Q 4 (B,x,t) 

Qs(B,x,t) 

Qe(B,x,t) 

Q 7 (B,x,t) 

Qs(B,x,t) 

k.pj(Xj't) ~ 

k gi (x,t) - 

k hl (x»t) = 

k f2 (x,t) = 
k g2 (x,t) = 
k h2 (x,t) = 
k f3 (x,t) = 
k g3 (x, t ) = 
k h3 (x,t) = 

Ri(x,t,a) 


_ B . 2B(3t-x)(t+x) 4Bxtf3(t+x) 2 -B 2 1 

” B z +(t + x) 2 [B 2 +(t+xj 2 J 2 [B 2 +(t+x) 2 ] d 

_ BrB 2 -(t-x) 2 l 

- [BMt-xFJ? • 

- B . 2B(t+x) 2 4Bxt[3(t+x) 2 -B 2 ] 

" B z +(t+x) 2 [B z +(t+x) z ] z [B^+lt+x) 2 ]" 5 

= (x-t)[(t-x) 2 -4B 2 ] 

L4B 2 +(t-x) 2 J 2 

_ (t-x)[(t+x) 2 -4B 2 1 4xt(t+x)r(t+x) 2 -l2B 2 l 

' [ "4 B 2 +(t+x7 z p [4B‘^t+x)’ z J 3 • 


(B.3) 


/ [Ri(x,t,a)+R 2 (x,t,a)-R 1 (H-x,H-t,a)-R 2 (H-X,H-t, a )]da, 

f° [R 1 (x,t,a)+R 2 (x,t,a)-R 1 (H-x,H-t,a)-R 2 (H-x,H-t,a)]cos(aB)dct, 

f g^y{Si(x,t,a)+S 2 (x,t,a)+S 1 (H-x,H-t,a)+S 2 (H-X,H-t,a)]sin(aB)da, 

f Q^y[Ri(x,t J a)+R 2 (x,t,a)-R 1 (H-X J H-t,a)-R 2 (H-X,H-t,a)]cos(aB}da, 

2 /^^j[Ri{x,t,a)+R 2 (x,t,a)-Ri(H-x,H-t,a)-R 2 (H-x > H-t,a)]cos 2 (aB)da, 

/ ■Q^y£s 1 (x,t,a)+S 2 (x,t,a)-S 1 (H-x,H-t,a)-S 2 (H-x,H-t,a)]sin(2aB)da, 

00 1 

f Q^jCTi(x,t,a)+T 2 (x,t,a)+T 1 (H-X,H-t,a)+T 2 (H-x,H-t > ct)]sin(aB)da, 

/ 0 ^y[Ti (x,t ,a)+T 2 (x,t,a)+T 1 (H-X ,H-t,a)+T 2 (H-x,H-t,a)]sin(2aB)da, 

2 /^yCU 1 (x,t,a)+U 2 (x,t,a)-U 1 (H-x,H-t,a)-U 2 (H-X,H-t,a)]sin 2 ( a B)da, 

(B.4) 


= [4a 2 H 2 -e (2-3at-aX+2a 2 Xt) + 2-3at-aX+2a 2 xt 

-2a 2 H 2 + 2a 2 tH + 2a 2 xH - 4aH]e' a ^ t+x) , 
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R 2 (x,t,a) = [(l-e' 2aH )(aX-at-2)-2aHCctX-2+2a 2 xt-2a 2 xH-3at+3aH)]e a(t " x) , 

(B.5) 

S 1 (x > t»a) = [ (4a 2 H 2 -e" 2a ^) ( 1 -aX+3at-2a 2 xt )+l -ctX+3at-2ct 2 xt+2a 2 H 2 
-2a 2 xH + 2aH - 2a 2 tH]e' a(t+X ^ , 

S 2 (x,t,cc) = [(l-e‘ 2aH )(aX-l-at)+2aH(aX-l+2a 2 xH-2a 2 xt-3aH+3at)]e a(t “ x) 

CB.6) 

Tl(x,t,a) = [(4a 2 H 2 -e“ 2aH )(-l+aX+at-2a 2 xt)+2a 2 H 2 -2a 2 xt-l+aX+at 
+2aH-2a 2 xH-2a 2 tH]e" a(t+x) , 

T 2 (x,t,a) = [(T-e~ 2aH )(at+l-ax)+2aH(aX-l+2a 2 xt-2a 2 xH-at+aH)]e a ^ X ^, 

(B.7) 


Mx.t^a) = [(4a 2 H 2 -e’ 2aH )(aX+2a 2 Xt-at)+aX+2a 2 Xt-at-2a 2 H 2 +2a 2 xH 
+ 2a 2 tH]e" a ( t+x ) , 

U 2 (x,t,a) = [0-e~ 2aH )(at-ax)-2aH(aX+2a 2 xH-2a 2 Xt-aH+at)]e ct ^ t X ^ , 

(B.8) 


k 14 (x) = . -^y [Vj (x,0) + V 2 (m,x,0) + V 2 (H-m,H-x,0)] , 

k 24 (x) = - [V 1 (x,B)+V 1 (x,-B)+V 2 (m > x > B)+V 2 (m,x,-B) 

+ V 2 (H-m,H-x,B) + V 2 {H-m,H-x,-B) + k pl (x)] , 
k 34 (x) = - ^ ■ < 1 fT y [Y 1 (x,B)-Y 1 (x,-B)+Y 2 (m,x,B)-Y 2 (m,x,-B) 

-Y 2 (H-m,H-x,B)+Y 2 (H-m,H-x,-B)+k p2 (x)] 


.. , (<+3)(B-n) + 2(x- m) 2 (B-n) 

V 1 (x,B) ~ - 2[.(x-m) z + (B-n) z J [(x-m) z + (B-n) z J z 


(B.9) 
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(B.10) 


v 2< ra > x ’ B > = ~ 2 [ ix^ ( MB-n)^J * 

- ■ 

[b'^FP ’ 


Y,(x,B) - - 


U+3)(x-m) . 2(x-m) 3 

2L(x-m)^+(B-n ) 2 J L(x-m) z+ (£ 




2 -(B-n) 2 ] 


4mx(x+m 


t{ 


2 -3(B-n) 2 ! 


x+m] z +[B-n 


k (x) = / [Xi(m,x,a) + X 2 (m,x,«) + X^H-m.H-x.a) 

P 1 o 

+ X 2 (H-m,H-x,a)]sin(an)sin(aB)da, 

GO 

k p2 (x) = / [Z 1 (m,x,a)+Z 2 (m,x,a)-Zi(H-m,H-x,a) 

P o 

-Z 2 (H-m,H-x,a)]sin(ctn)sin(aB)da, 
Xx(m,x,a) = { [K-2am)(2ax-3)-l](e ^-4a 2 H 2 -1) 

+ 2aH(2am-2aH+2aX-3-K:)}e" a ^ m+x) , 

X 2 (m,x,a) = {(1-e (2am-2ax+ic+3) 

2aH[2aX ( K+2am-2aH ) -3ic+6aH-6ani- 1 ] }e a (m “ X J 


Zi(m,x,o) = { (l-e“ 2aH +4a 2 H 2 )[(ic-2am)(2ax-l )-l] 

+ 2aH(2aH-2am-2ax+2+K) }e a ^ m x ^ , 

Z 2 (m,x,a) = {(1-e (2ax-2am-K-l ) - 2aH[2ax(ic+2am-2aH) 
- K - 1 + 2aH - 2am]}e a(n, " X ^ . 


(B.ll) 


(B.12) 


(B.13) 


(B.14) 
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Figure 1 . 


The basic crack geometry. 



Figure 2. Stress intensity factors and probable crack propagation 

angle in an infinite strip containing two edge cracks under 
uniform tension, d=0.2H. 
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Figure 4. 


Stress intensity factors and probable crack propagation 
angle in an infinite strip with two edge cracks under 
bending, d=0.2H. 
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Figure 6. Stress intensity factors in an infinite strip containing 
three edge cracks under uniform tension, d=b=0.2H. 
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Stress intensity factors in an infinite strip containing 
three edge cracks under bending, d=b=0.2H, o m =6M/H 2 . 


Figure 8. 
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Figure 11. Stress intensity factors in an infinite strip containing 
edge cracks and subjected to three point bending, a m = 
6M/H 2 = 24P/H. 
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Figure 12. Stress intensity factor in a square plate containing a 
wedge-loaded central crack. 
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Figure 13. Stress intensity factor in a uniformly loaded rectangular 
plate naving an edge crack. 



Figure 14. Stress intensity factor in a square plate with an edge 
crack under wedge loading. 
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